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SUMMARY 


The differential equations of the bending of sandwich plates were 
integrated to obtain the deflections when the four edges of the plate 
are simply supported and the loading consists of either a uniformly dis- 
tributed transverse load or a concentrated load applied at the center 
of the panel. The deflection patterns are shown in diagrams and the 
maximum deflection of the plate is presented in a number of graphs. 


INTRODUCTION 


A sandwich plate is a composite plate consisting of two thin faces 
and a thick core. In airplane construction the faces are usually com^ 
posed of aluminum alloy, and the core is congaosed of some lightweight 
material such as an expanded plastic or balsa wood. In the latter case, 
the fibers of the wood are usually arranged perpendicular to the plane 
of the plate. Since the modulus of elasticity of the core is of the 
order of magnitude of one-thousandth that of the faces, the normal 
stresses in the core are of little importance in resisting bending 
moments, although the usual ratio of face thickness to core thickness 
lies between one-tenth and one-hundredth. The core performs a task in 
transmitting shear forces and undergoes considerable shearing deforma- 
tions because its modulus of shear is low, and therefore shearing defor- 
mations cannot be disregarded in the analysis of sandwich plates. 

Differential equations have been derived for rectangular sandwich 
plates subjected to transverse and edgewise loading (reference 1) . In 
the present report the differential equations are integrated for the 
case when all four edges of the plate are simply supported, and when 
the load is either concentrated at the center or unifonnly distributed 
over the entire plate. The maxi, mum deflection depends upon the thick- 
ness ratio r = c/t and upon a nondimensional parameter R, Numerical 
values of the deflections of a square sandwich plate were calculated for 
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a great member of values of these two parameters, and the results of the 
computations are presented in the form of diagrams. 

The calbulations presented here were carried out under the sponsor- 
ship and with the financial assistance of the National Advisory Committee 
for Aeronautics. The authors are indebted to Doctors N. J. Hoff and 
V. L. Salenio for their advice in the course of the calculations and for 
their help in preparing the final report, and to Mr. George Booth for 
his work in carrying out the calculations and in checking the analysis. 


SYMBOLS 


B 

B 


c 

d 


c 


D 

Df 

Do 

E 

F 

Do 

k 


Lx 


deflection factor for sandwich plate under concentrated load 

deflection factor for sandwich plate xmder uniformly distrib- 
uted load 

thickness of core, inches 

bending rigidity of thin plate, pound-inches squared per inch 

bending rigidity of two independent faces, potind-inches 
squared per inch 

bending rigidity of sandwich plate, pound-inches squared 
per inch 

Young's modulus of face, psi 
form factor for sandwich plate 
shear modulus of core, psi 
side ratio 

side length of sandwich plate in x-direction, inches 

side length of sandwich plate in y-direction, inches 
concentrated load, pounds 

con 5 )ressive end load in x-direction, pounds per inch 
compressive end load in y-direction, pounds per inch 
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q distributed load, psi 

intensity of uniformly distributed load 

r = c/t 

R stiffness factor for sandwich plate 

t thickness of face, inches 

u displacement in x-direction, inches 

V displacement in y-direction, inches 

w deflection in z-direction, inches 

x,y rectangular coordinates in plane of faces, inches 

z rectangular- coordinate perpendicular to plane of faces, inches 

[I Poisson's ratio 

p 

A Laplace operator 

DERIVATION OF EXPRESSIONS FOR DEFLECTIONS 

The differential equations for the deflections of a sandwich plate 
have been derived previously (reference l), and the problem has been 
defined by means of the following three partial differential equations 

® 0 1^^ + (1 - + (1 + " 2GqCu - 2GcC ° ^ ^ "x = ° 

+ (1 + " 2GqCv - 2GqC 'V " ° 

Df - GcC — ^ (u^ + Vy) - G^c A^w - q + P^w^ + P^Wy^ = 0 (Ic) 

together with the following boundary conditions: 

u^ + pvy = 0 when x = 0, 


(2a) 
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II 

when 

7 = 0, 

V 

(2b) 

0 

when 

7 = 0, 

L 

7 

(2c) 

0 

when 

II 

O 

Lx 

(2d) 

= 0 

when 

o 

II 


(2e) 

= 0 

when 

7 = 0, 

V 

(2f) 

0 

when 

X = 0, 

and when y = 0^ Ly 

(2g) 


The symbols denoting bending rigidities were defined in the following 
manner: 


Dq = Et(c + t)^y^[2(l - \s?)_ 


(3a) 


Df = Et3/[6(l - (3b) 

where D^. is the bending rigidity per inch of the faces about their own 
centroidal axes, calculated for the two faces. Also, is the bending 

rigidity of 1-inch width of the sandwich panel calculated about the 
centroidal axis of the sandwich, when the contribution of the core, as 
well as that represented by D^, is neglected. 

The other symbols, as well as the sign convention, are shown in 
figure 1. It should be mentioned, that the sign of the last term in. 
equations (la) and (lb) and that of the second term in equation (ic) are 
opposed, to those in reference 1 because of the different choice of the 
coordinate systems. A single sixth-order differential equation can be 
obtained, from the above three differential equations by elimination of u 
and V (reference l) : 


A^w 


- (°o Df)(G^,c/Do)A^w = - (g^c/Do)^ (q - P^w^ - PyWyy) 


(U) 


^6 » {d^/dx.^) + 3{d^/d^dy'^) + 3(ay3x^3y^) + {d^/dy^) 0;a) 


where 
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This problem will be solved by the method of Fourier series 
expansion. The deflection w can be represented by a double Fourier 
sine series: 


w 


OO 00 

m n 


Wmn sm 


^mra/Lx) sin (nny/Ly) 


( 5 ) 


which satisfies all the boundary conditions (2e) to (2g). Assungjtion of 
the other two deflections in the form 


and 


OO 00 

u = ^ ^ u^jjn cos (mnx/Lx) sin (nny/Ly) 


00 OO 

V = ^ ^ Vjjjn sin (miuc/Lx) cos (nny/Iy) 


m n 


(6a) 


(6b) 


satisfies all the remaining boundary conditions (equations (2)). Equa- 
tion (2d) implies Yj = 0 on x = 0 and x = L^, and equation (2a) 

then reduces to u^^ =0 on the edges x = 0 and x = Likewise^ 

equation (2c) implies u^ = 0 on y = 0 and y = ly, and equation (2b) 
then reduces to vy = 0 on the edges y = 0 and y = ly. If q is 

expressed by means of the Fourier sine series 


00 OO 

^ - Z Z sin ^miix/L^j.) sin /nny/ly) (?) 

m n ' 

where the q can be calculated from the well-known formula 

^mn 


- (VVy) 



( 8 ) 


then the Fourier coefficients and Wnm can be determined 

from the algebraic equations obtained by substituting u, v, w, and 
q into equations (la), (lb), and (ic).' The knowledge of the coeffi- 
cients Ujjj^ and is not needed when the transverse deflection of 

the face plate is sought. Therefore, it was found more convenient to 
solve the problem by means of the sixth-order differential equation 
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(equation (ii.)). Substitution of q from equation (7) and w from 
equation (5) into equation (i;) yields: 



In the present solution of the problem^ the compressive end 
loads and will be assumed to be zero. With the introduction 

of the following notations 


k ■ (10a) 

r = -c/t (10b) 

D = Ety[l2(l - ix^)] (lOc) 

R = GgtLyyn^D (lOd) 


the deflection coefficients (equation (9)) can be written in the fol- 
lowing nondi mensional form: 


w 


mn 


mn 




( 11 ) 
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where 


Cjnn = (i/2k^) |f6(l + r)Vr^(m2 + (lla) 


and 


Kjju^° = [6(1 + r)yrli](m2 + + [3(1 '+ r)2 + l](m2 + (lib) 


The maxiimm deflection occurs at the center of plate. Substitution 
of X ■ Lx/2 and y * Ly.y'2 into equation (5) yields 


w 


max 




w (_l)l+C(ni+n)/2l 

mn' 


(12) 


where w^^^^ is given by equation (U). 


UNIFORMLY DISTRIBUTED LOAD 


When the load is uniformly distributed over the plate, the coef- 
ficients of the Fourier series given by equation (?) can be calculated 
as 


- (K/Vy) 



(mnx/L^) sin ^nny/Ly.) dx dy 



(13) 


where q^ is the intensity of the uniformly distributed load. 

Substitution of these coefficients into equation (ll) and intro- 
duction of the notation 


B 


d 



(ih) 
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yield the following ejqjression for the nondimensional deflection 
parameter: 


00 CO . 

H^d) = 21 Z_ i^mn/^) sin sin ^niiy/Ly) 

m=l,3,‘" n=l,3, ••• ' 


(15) 


The maximum deflection parameter vfmayl^d can he written as 


(w^/B,) = (l6/n2) ^ 

^ m=1^3,*** n=1^3, ••• ^ 

The maYiimnn deflection can be calculated from the rapidly convergent 
series (16) • 


to: 


If c or approaches zero, the ma^ylTmim deflection is reduced 


w, 


max 


= K/A) ^ £ 

' m=l,3, ••• n=l,3,”' 


(- 1 ) 


l+[(m+n)/22| 


mn 


]m/l^)2 + (n/Ly)2]' 


(17) 


This is the formula for the maxi mum deflection of a s imp ly supported 
thin plate under the uniformly distributed load qo/^-» given in 

reference 2. In the limiting case of no shear deformation in the core 
(infinite Gg), the maximum deflection is given by equation (17), with 
D replaced by DF. The form factor for the sandwich plate is 
F = 1 + 3(1 + c/t). The quantity DF represents the total moment of 
inertia of the two faces with respect to the center plane of the sand- 
wich plate . 

Deflections of sandwich plates under uniformly distributed load are 
shown in figures 2 to ij.. 
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CONCENTRATED LOAD 


In the case of a concentrated' load P acting at the center of" the 
sandwich plate, the Fourier coefficients in equation (7) are 

obtained by a limiting process. Let p be the intensity of the load 
uniformly distributed over a small square with sides a parallel to 
those of the plate, and with its center at the center of the plate. 

Set P = pa^. Then according to equation (8), 


%n 





(p/a^) sin ^mnx/Lj.^ sin dy 


If this expression is integrated and the quantity a is allowed 
to approach zero while P = pa^ is kept constant. 


'‘mn 


when m and n 

when m. and n are even. 

With the notations used before, and with 


are odd, and qjjjj^ = 0 



( 18 ) 


the deflection can be expressed in the following nondimensional form: 

* 00 oQ X’^l r 771+11 ) / 2 I 

(w/Bc) =idc 2T ^ Cmn("l) sin (mjtx/Lx) sin (nity/Iy) 

m=l,3,**’ n=l,3j**' 


(19) 


The maxj.mum deflection is then given by 


00 00 

(wmax/Bc) Idc ^ 2_ 

^ ' m=l,3,--- n=l,3,-‘ 


mn 


(20) 


where Cjjjj^ is defined by equation (Ua) . 
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In equation (20) if either c or G^, approaches zero, the maxi- 
mum deflection reduces to that of two thin plates with no core (refer- 
ence 2), namely: 


(«/Bc) - (seA**) 


S ' 




( 21 ) 


It is shown in equation (a 8) of the appendix that the value of the 
double series for a square plate (k = l) is 


m=l,3,'” n=l,3,-'- (m2- + n2)‘ 


= 0 . 282^1 


(21a) 


Therefore the maximum deflection of the two tMn square plates with 
no core is 


("max/®c) = 5.800U2 x lo"^ (2ib) 

= a 

For a square sandwich plate the expression for the maximum deflec- 
tion given by equation (20) may be transformed as follows: 


00 00 

(wjnax/^c) = ct + 1; 'y y 

' m=l,3, ••• n=l^*** 




^ 2: 0 .- 

m=l, 3 , n=l,3, ••• 


m= 


=1,3, h=l,3,“* (m^ + n^k^) 


( 22 ) 


0 
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(W®c) =5.800i;2 X 10-3 _ 



nl,3,’-‘ + + + {xB./2) + [^rR/6(l + r)^ 

(23) 


R 


For a square sandwich plate with an infinitely rigid core (i.e., 
■> ”) the ma-xi Timm deflection, obtained from equation (20), is: 


(w^/Bc) - y-'') t: ti 

m=l,3, • • • n=l,3. 


n=l,3, *’* (ni^ + n^) 3(1 + p)2 + 


= (aA) 


where 


F = 3(1 + r)^ + 1 


( 21 ;) 

(2l;a) 


is the form factor of the sandwich plate. 

Another expression for the maxj.mum deflection of a square sandwich 
plate may be obtained from equation (20); 


('^max/^c) = ^ 

m=l,3,*'* n=l,3, ••• 


= 5.8001;2 X (io“3/f) + (2/n^)[(F - l)A] 




1 



•n2)(. 


2 2 
m*^ + n^ + 


(rR/2) + 


|rR/6(l + r) 



( 25 ) 
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Equations (23) and (25) are alternate expressions for the deflec- 
tion at the center of the sandwich plate. Since the series involved 
are sums of positive constants, the following inequa l ity holds: 


“ > ^max/Bc > (^6) 

This means that the deflection of a sandwich plate is bounded by two 
limiting cases: It is smaller than that of a sandwich plate with a core 

having no shearing resistance and is greater than the deflection of a 
sandwich plate with a core of infinite sheaidng resistance. 

Equation (23) can be used for values of R < 1 and r < 50 since 
the convergence of the series is sufficiently rapid for these values. 
This series was used to obtain portions of figure 5. ^or large values 
of R the convergence of both series in equations (20) and (23) was 
too slow to permit direct use. Because of its sirgjler algebraic form, 
the series of equation (25) was evaliiated as shown in the appendix. The 
series (A12), in conjunction with- other modifications given in the 
appendix, was used to obtain figure 6. 

PRESENTATION OF RESULTS 


The equations derived in the preceding section have been used to 
compute the deflections of various square sandwich plates subjected to 
unifomOy distributed loads or to a concentrated normal load acting at 
the center of the plate. 

The series solution equation (l5) was used -bo obtain the maximum 
deflection of a square sandwich plate supporting a uniformly distributed 
load. The deflections for a square sandwich plate at y = Ly.^2 and for 

different values of x/L^j- have been obtained for R = 800 and dif- 
ferent values of r. These are plotted in figure 2, which shows the 
deformation pattern for the sandwich plate. The maxi mum deflection for 
a sandwich plate with prescribed values of r and R is shown in fig- 
ure 3. In figure ij., the maximum deflection for values of R from 0 

to 1 is plotted. As R approaches zero, the deflections approach the 
value for -a sandwich plate without shearing deformation, that is, 

= 2.029 X lO”^ inches. Hence the curves of figure 1; show the 
effect of the core upon the maxi mum deflection of a sandwich plate. 

For a square sandwich plate with a concentrated load at the center, 

the maximum deflection for small values of R is plotted in figure 5* 

In figure 6 the maximum deflection for large values of R is plotted. 
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In both the uniformly distributed and the concentrated-load cases^ 
as R approaches infinity the deflections approach the asymptotic 
values of a sandwich plate without shear deformation. The limiting 
values for a square plate are given in figures 3 and 6 for these two 
cases, respectively. 

With the aid of figures 3 and 6, the maximum deflections of a 
square sandwich plate can be calculated. The value of may be 

read off from the curves of figure 3 for appropidate values of r and 
for values of R less than lOCXD. If R exceeds 1000, figure 3 serves 
to obtain only an approximation to the value of . In such cases 

only three terms of equation (l5) are needed to obtain v^max/^d s-ccu- 
rately, and in most applications the first term of equation (l5) would 
be adequate. For values of R less than 1000 the graphs of figure 3 
are particularly valuable since in this range many more terms of equa- 
tion (l5) would be required to calculate accurately. 

COMPARISON WITH EXPERIMENT 


The calculated values of the maxi.mum deflection under uniform 
loading have been compared with the test results obtained by the Forest 
Products Laboratoiy (reference 3f table l). The results calculated from 
the present report for a uniformly distributed load are in good agree- 
ment with these test results, although of slightly smaller value in 23 
of 30 cases. This behavior is to be expected since the theoretical 
boundary conditions correspond to a greater degree of constraint along 
the edges than do the experimental boundary conditions. Theoretically, 
both faces must be simply supported; experimentally, only one face was 
simply supported on knife edges. 


This comparison is shown in tables I, II, and III. The average 
absolute errors are k-7, h-9, and 6.5 percent, respectively. In l8 of 
30 cases the absolute error is less than 5 percent; in one-third of the 
cases the absolute error is less than 3 percent. The four largest per- 
centage errors occurred for the case in which the 'core was composed of 
corrugated paper honeycomb. 


NDMERICAL EXAMPLE 


The following examples show how to find the maximum deflection of 
a simply supported square sandwich plate by means of figures 3 and 6. 



Given: 


Ljj. = ly = 10 inches 

c = 0.65 inch 
t = 0.025 inch 

E = 10.5 X 10^ psi 
H = 0.3 
Gg = 6000 psi 

First calculate the following parameters from equations (lOb), 
(10c), and (lOd): 

r = c/t = 26 

D = EtYl2(l - = 15.0 

R = GctL^yn^D = 101 


Uniformly Distributed Load 
For the uniformly distributed load, from figure 3, 

("mx/Bd) " 3.9 X 10-6 

For a uniformly distributed load of intensity = 1 psi the 
value of according to equation (lii) is 

Bd = qoBxV® = 

Hence the ma-xn irmni deflection is 


w. 


33iax 


= 0.00260 inch 
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Concentrated Load at Center of Plate 

Consider the concentrated load equal to the total uniformly dis- 
tributed load in the section above: 

P = = 100 

Be = 

From figure 6 

("max/Bo) ' 28 x 
Hence the maximum deflection is 


Bd 

667 


"^max = 


CONCLUDING EEMEKS 


The deflections of a square sandwich panel were calculated for the 
case where all four edges are simply supported and the external loading 
is either a uniformly distributed transverse load or a concentrated 
transverse load applied at the center of the plate. Diagrams are pre- 
sented which show the pattern of deflections for these cases. The main 
results of the calculations are the graphs from which the mayimum 
deflection of the sandwich plate can be read as a function of two non- 
dimensional parameters. 

The faces were assumed to be isotropic thin plates. The distance 
between the faces was assumed to remain -unchanged during the defor- 
mations, and the contribution of the core to the bending rigidity of the 
entire plate was neglected. Theoretical as well as experimental con- 
siderations have shown that these assumptions give satisfactory results 
for sandwich plates of the types in use today in the United States. 

Satisfactory agreement with the experimental results of the Forest 
Products Laboratory was obtained. 
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APPENDIX 


ON EVALUATION OF DOUBLE ' SERIES S(a,b) 
Introduction 


The eqxiatdjOn for the double series S(a^b) is: 

S(a,b) = ^ + a^) (m^ + + b^^ (Al) 

m=I,3,--- n=l757*” 

The value of the series S(a,0) was required to a high degree of 
accuracy in order to evaluate the deflection equation (25) developed in 
the bo(^ of this report. The main set of values considered for a^ 
ranged from 200 to lj.0,000 (b = O) . In addition, the special value a = 0, 
b = 0 had to be considered. This latter case S(0, O) was treated 
separately. 

Direct summation of a finite number of terms of equation (Al) does 
not yield the required six-figure accuracy in a reasonable time, par- 
ticularly for the larger values of a (b = O) . Methods were developed 
to determine S(a,b) quickly to any degree of accuracy. 


EvalTiation of S(0,0) 

For the case a = 0, b = 0, equation (Al) becomes 


s(o,o) = 51 51 

m.=l,3, ••• n=l,3,-*' 


(A2) 


The series (A2) can be summed explicitly by rows by using the 
partial-fraction expansion of the hyperbolic tangent (reference U, 
p. 362); 



2 + ^2 


) 


(tx/U^) tanh (n^l/2) 


(A3) ■ 
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Differentiating both sides of equation (A3) with respect to t 3 and 
replacing i9 by ‘ m in the result lead to 


n=l,3. 


l/(m2 t n2)' 



( nm/2 ) sech 


2 



(All) 


The result (Aij.) permits the summation with respect to n to be 
carried out in equation (A2); S(0,0) may now be written as the single 
series 


S(0,0) = (n/8) 


m=l. 


{l/m^ 




tanh (nm/2) - (;tm/2) sech'^ 



It would be difficult to treat equation (A^) explicitly in the 
form shown, but the evaluation can be simplified if the dependence of 
the hyperbolic functions upon m is considered. 

For example, with m = 7, 


tanh (7n/2) = 1 - 6.7k x 10"^ 
sech^(7n/2) = 1.35 x IQ-^ 


1 


(A6) 


If m > 7^ then 1 - tanh (nm/2) and sech^(nm/2) approach zero 
even more rapidly. For the accuracy required this shows that 
tanh (nm/2) may be set equal to unity and sech^(nm/2) may be set 
equal to zero for m ^ 7. This singslifies the form of equation (a5).* 

If the series 


S' 



(A7) 


can be evaluated, this value can be used to determine equation (a 5) to 
any degree of accuracy, for only the first three terms of equation (A7), 
multiplied by n/8, need be subtracted. To this value must be added 
the first three terms of equation (a 5). The result will be S(0,0) to 
six significant figures. 
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The series (A?) must be calculated directly. It is possible to 
transform series (A?) into a rapidly convergent series (reference 5, 
p. 272), from which S' may be determined to be 

S* = 1.051 799 8 ••• 

The value of S(0,0) may now be detennined; 

S(0,0) = 0.282 51 • • • (A8) 


Evaluation of S(a,0)j (b = O) 
Equation (Al) may be written in the form 


S(a,0) = (l/a^) J [V(“^ “ V^( 

m=l 737 • • • ^=1757 ’ * • ^ 


2 2 2 


[ 


(A9) 


Application of series (A2) to series (A9) reduces the double 
series (A9) to the single series 


S(a,0) = (n/ka.^) 


m=l,3, ••• 


(l/m) tanh (nm/2) - 


(l/\|m2 + 


tanh 



(■AlO) 


The same approximations (A6) which s impli f i ad series (a 5) into (A7) 
may be -used to simplify series (AlO)j the tsrpical results (A6) show that 
the hyperbolic tangents may be replaced by unity after the third term. 
Since the smallest nonzero value of a is approximately li;, the second 
hyperbolic tangent of series (AlO) can be set equal to unity in the 
first term. The series (equivalent to series (A7) in the case b = O) 
which must now be considered is 


S*(a,0) 


y~ (l/m) - ^1 y^\|m^ + a^^ 


m=l,3. 


(All) 
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If S*(a,0) of series (All) were known, then S(a,0) of series (AlO) 
could be determined quickly, as before, by taking into account sepa- 
rately the first three terms of series (AlO). The series (All) is not 
suitable for accurate and rapid numerical calculation. 

The value of S*(a,0) of series (All) may be determined to a high 
degree of accuracy by means of the following result, which will be 
proved. 


S*(a,0) = I y + I loge a + JL ("l)^"^^ Ko(nna) (a12) 

n=l727 • • • 

where y = Euler's constant (y = 0.^77216 • • •) and Kq is the modified 
Bessel function of the second kind of order zero. Because of the prop- 
erties of Kq and the range of val,ues of a, it will be shown that 
the infinite series in equation (A12) can always be neglected, its 

largest value being of the order of magnitude 10“^^. Then series (A12) 
can be used to evaluate S^(a,0) to a high degree of accuracy and 
S(a,0) can then be determined quickly with the known value of S*(a, O). 

The result (A12) will be established by converting equation (A3J.) 
into an infinite integral and evaluating the result. A check upon the 
work is possible by means of a contour integration; this will be dis- 
cussed briefly later. 

To convert equation (All) into an integral representation, the 
following result from the theoiy of Bessel functions is needed (refer- 
ence 6, p. 6^) 



Jo (at) dt 



(a13) 


where is the Bessel function of the first kind of order zero; 

another result needed is 


e-mt dt = 1/] 


m 


(AlU) 


0 
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obtained by setting a = 0 in equation (A13), or by direct means. 
Combination of equations (A13) and (Alii) yields 


(1/m) - 




|l - Jo(at^ e""^^ dt 


(Al^) 


Substitution of equation (Al^) into equation (All) yields 


S*(a,0) = 


ni“l ^ 3 ^ • 


Cl/m) -[x/\[ 


+ a^ 


loo / 00 


0 \m=l,3,-- 


,-^(l - j^(at)) 


dt 


(A16) 


Now 


m=l,3. 


y e = e ^ ^1 e~^^ + e + • • •) 


= (e “ e (summing the geometric series) 

= ly/(e'^ - e-"*^) 

= 1/(2 sinh t) (A17) 


Insertion of equation (Al?) into equation (Al6) yields 


S*(a,0) 


1 P” [1 - Jp(at)] 

2 J Q (sinh t) 


(A18) 


Equation (A18) is the required integral representation of S*(a,0). 
This integral can be sin^jlified by using the result (reference 7, 
p. 136): 


^ . 1 . 2t 

sinh t t 


oo 

n‘1,2,- 


y(t2 * „2„2) 




(A19) 
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Insertion of equation (A19) into equation (Al8) yields 


S*(a,0) = i fl - (V-t) - 2t f (-l)^‘^y(t^ + n^n^) 

2 Jo L n=l^--* ' 


dt 


(A20) 


One of the typical integrands to be considered in equation (A20) is 
of the form tJQ(at)y/(t^ + but 




2 2 "\ 
n 


dt = RgCnna) 


(A21) 


is a known result (reference 6, p. ?8), where Kq is the modified 

Bessel function of the second kind of order zero. This fTinction is 
tabulated, for example, in Watson's "Bessel Functions" from 0 to l6 at 
intervals of 0.02 (reference 8, table II, pp. 698-713). The asymptotic 
formula for K^(x) is (reference 6, p. 55): 



i2/1!(8x2 + Ql^ X 3^)/21(8x)^ - 



(A22) 


The usual series representation (reference 6, p. 22) may be used for 
small, values of x. 


If the smallest value of a^ is 200, then a is approximately lij,. 

For the case n = 1, = e“^ = 10’"^^ (approximately) . This is 

the largest esqjonential term that can occur and shows that all contri- 
butions of the form of equation (A2l) will be entirely without influence 
if numerical work to six or seven places is required. However, if a 
were small (a = 1, or less) it would still be a simple matter to 
include terms such as those in equation (A21) . With the terms (A2l) 
considered negligible, equation (A20) may be rewritten as 


S*(a,0) = (1/2) <|(l/sinh t) 



dt 


(A23) 
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The integral in equation (A23) is far simpler than that in equa- 
tion (A18) since the term sinh t in equation (A23) has now been sepa- 
rated from IgCat). The integrals in equations (Al8) and (A23) converge 

at t = 0, and both infinite integrals exist. The integral (A23) can 
be treated as two integrals, one over the range 0 ^ t ^ l/a and another 
over the range l/a ^ t < ». Let these integrals be S^* and S2*, 
respectively; then 


S2* = (1/2) |^(l/sinh t) - jjQ(at)/^dt 

= (V2) ^ Ql/sinh t) - (1/t)] dt + (1/2) |^[l - Jo(at^/ 


tj dt 
(A21;) 


The term l/t was added and subtracted in equation (A2h) so that both 
integrals still converge at t = 0. After an integration by parts, the 
second integral of equation (A21;) becomes (using dJQ(x)/dx = -j^(x)) 



J^(at) logg (at) 


t=l/a 


t=0 



J]_(at) logg (at) dt 

(A25) 


The first term of equation (A23) is zero at both limits. The first inte- 
gral of equation (A2I4.) may be written as 


1 im 


Pl/a 

(1/2) J j^(l/sinh t) - (l/t^ dt = 


lim (1/2) log Itanh (t/2)/tl^ 

e 0 ■- -'■t=6 


which is equal to 


(1/2) logg [(tanh l/2a)/(l/a^ - (I/2) logg (l/2) 
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since the limit of (tanh e )/e is unity as e approaches zero. Then, 


= (1/2) logg tanh (l/2a) + (l/2) logg a + 


'1/a 

(1/2) logg 2 - (a/2) j J 2 _(at) log-g (at) dt (A26) 


The remaining integral is 




((l/sinh t) - 




(l/sinh t) dt “ 



[jo(a't)/£j 


dt 


(A27) 


Since both integrals of equation (A2?) converge separately • at the 
lower limits, the integrals may be considered separately. Integration 
of the second integral of equation (A2?) ty parts leads to the result 


(-l/2)Jg(at) logg a-^ 


t=oo 

t=l/a 



J]_(at) logg (at) dt 


(A28) 


The first term of equation (A28) is zero at both limits. The first 
integral of eqxiation (A27) is 


■b=oo 


_(l/2) logg tanh (t/2^ _ 


t=l/a 


(-1/2) loge tanh (l/2a) 


Addition of S^* and S 2 * yields 


S*(a,0) = Si* + S 2 * 

poo 

= i logg (2a) - (a/2 ) J Ji(at) logg (at) dt (A29) 
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The terms containing tanh l/2a cancel and the two integrals can 
be combined into one integral. Now (reference 6j p. 76), 



J]^(at) logg (at) dt = / Ji(x) log^ x dx 


= -7 + logg 2 (setting x = at) (A30) 

where 7 = Euler's constant = 0.^77 216 • • Inseirbion of equa- 
tion (A30) into equation (A29) and simplification yield 


S*(a,0) = (1/2)7 (1/2) logg a (A31) 

Terms of the type in equation (A2l) have been neglected in equation (A3l) . 
The exact result (A12) is valid for all a > 0 and is repeated here for 
convenience: 


S^a,0) = Y1 

m=l,3, • * * L_ 


(1/m) - (l/{ 


2 2 

m^ + ar 


= (1/2)7 + (1/2) log a + Y1 Kg(nna) (A32) 

n=l,2 , ••• 


General Case S(a,b) for a / 0 and b / 0 

The use of equation (A32) permits S*(b,0) - S*(a,0) to be written 
as (a > Oj b > 0) : 


y (-1)"'"^ rKo(nnb) - K (nna)1 

n=l^,-*- -■ 


= (1/2) loge (b/a) + 


(A33) 
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The series (A33) would be needed in the evaluation of the general 
series S(a,b), which is a natural generalization of equation (25). It 
is now possible to evaluate equation (A33) in precisel'-’’ the same way as 
the series (Al) was evaluated for the case of a 0, b =■ 0. The case 
of a = Oj b = 0 has already been treated. Only the case of a = b 0 

must be considered to con^jlete the discussion of equation (A33). Equa- 
tion (A33) can be written as 


S(a,b) 








(A3li) 


The auxiliary series for equation (A3U), which corresponds to eqiia- 
tion (All) in relation to equation (A9), is 


S^a,b) = l/(b^ - a^) ^1 /\|m^ + a^ j - ^1 j \Jm^ + b^ j (A35) 


Application of equation (A33) yields 


(b + a) 1 b - a 


(- 1 ) 


n+1 


n=I737 • • • 


KQ(nnb) -KQ(nna) 


b-a 


(A36) 


It is possible to let b approach a in equation (A36) to obtain 


S (a,a) 


00 



(-1)^ K^(nna) 


(A37) 


since the expressions of equation (A36) define the derivatives of (logg a) 
and of Ko(nna), respectively, with dKQ(x)/dx = - K-j^(x). The func- 
tion K]_(x) is the modified Bessel function of the second kind of order 
one and is also tabulated in Watson' s "Bessel Functions" (reference 8, 
table II, pp, 698 - 713 ). Thus equations (A37) and (A33) can in general be 
used to evaluate equation (Al) for all values of (a,b) which are not 
both zero. 
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The result (A33) can be checked by a contour integration over the 
contour shown in figure 7j the details will be omitted. The proper • 
integral to consider is 


(lA) 




(az) - 




dz 


The contribution over the large semicircle approaches zero as R — > “ 
while along the real axis the result is 



which is the integral representation of equation (A33). The residue 
at z =■ 0 can be shown to be (l/2) logg (b/a) while the poles at 
z = nni (n=lj2, — ) contribute residues which are the terms of the 
infinite series (A33). This method shows why the original integral (Al8) 
could not be evaluated -in this manner, for it is not possible to let a 
or b approach zero in equation (A33). The case a = 0 must be treated 
separately. 

The method developed here for the evaluation of equation (Al) was 
not that used at the start of the calculations but was completed later 
and then used for all calculations. The method originally used may be 
s umma rized briefly as follows: Since equation (A3) permits an explicit 

s umma tion of equation (Al) by rows, a definite number of rows (actually, 1 $) 
was summed quickly by this method. Ey symmetry that is also the sum of 
the first l5 columns. To illustrate the procedure graphically (see 
fig. 8) let the entire quadrant represent S(a,0). Summation of the rows 
is represented by the shaded area 125678, and summation of the same 'num- 
ber of columns, by 123U58. The terms 1258 are counted twice and must be 
found separateHy. Thus, the sum of the terms contained in the L-shaped 
region 123li5678 is known; if ttiis sum is subtracted from the series 
itself, the remaining terms correspond to the unshaded area U56“>. This 
remainder of the series was approximated by the integral 



dx dy 

+ y^) (x^ + y^ 



where the lower limits u were chosen appropriately, depending upon the 
number of columns taken initially. This integral was used to determine 
upper and lower bounds for the remainder of the series. Much of the 
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work was directed toward obtaining a useful value of the integral for 
large values of a. In this way, bounds on S(a,0) were congjuted. This 
established an upper limit to the percentage error in the results. 

One such calculation led to the values (a^ = lj.0,000): 

Upper bound: 56.575 ^ 10“^ 

Lower bound: 55*912 x 10”^ 

The method using equation (A12) yields the value 56.0551 ^ 10”^ 
which lies between the calculated values of the upper bound and lower 
bound. 


Generalization of Series S(a,b) 

The series 


S(a,b,c) = y y * c^n^ + a^) + c^n^ + b^) (A38) 

m=l,3j“* n=l,3, ••• 

is the generalization of series (Al) to the plate of rectangular shape; 
c = 1 corresponds to the case of a square. Equation (A38) may be 
written in the form 

S(a,b,c) = l/(b^ - a^) ) ^ ^[l/(m^ + c^n^ + a^)"] - 

m=l,3,“- n=l,3,--- I*- 



(A39) 


The summation with respect to m may be carried out as before; the use 
of equation (A3) yields 


S(a,b,c) 


2 

U(b2 - a2) n-iX--- 


tanh (ti/2) \ 


c^n^ +. 


\ c^n^ + a^ 


tanh (n/2) \ c^n^ + b^ 


nI 


c^n^ + b^ 


(AiiO) 
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In this case, the auxiliary series is 


S*(a,b,c) = - 
c 


^ (l/\/n2 * * bVc^) 


(Ala) 


which can be evaluated by equation (A33). The case a = b can be 
treated as before by the use of equation ( A37 ) • Thus the general 
series (A38) can also be evaluated by the methods already developed. 


Pol 3 rtechnic Institute of Brooklyn 
Brooklyn, N. T., June 19, 1951 
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TABLE I 

COMPARISON OF TESTED AND CALCDLATED DEFLECTIONS OF A SIMPLY 

SUPPORTED SQUARE SANDWICH ELATE UNDER UNIFORM LOAD 

[^nmi num faces with E = 10 x 10^ psi and balsa core 
with G(j = 12,670 psi (reference 3)j average 
absolute error, ij..? percent] 


Plate 

• Thickness 
of faces 
(In.) 

Thickness 
of core 
(in.) 

Span 

(in.) 

Unit load 
on plate 
(psi) 

Test 

deflection 
(reference 3) 
(in.) 

Calculated 

deflection 

(in.) 

Percent 

error 

AB-U 

0.032 

0.?lli 

14i.20 

0,31;08 

0.168 

0.171 

1.8 

AB -3 

.032 



.390 

ijli .20 

.7211 

.359 

.377 

3.0 

AB-i; 

0.032 

0.301 

38. Oil 

0.7211 

0.123 

O.I 3 I; 

8.9 

AB -3 

.032 

.381 

38 . 0 ii 

.^90 

.189 

.198 

1;.8 

AB- 1 ; 

0.032 

0.307 

32.0ii 

0.7211 

0.072 

0.069 

-h>2 

AB -3 

.032 

.383 

32.01; 

.7211 

.112 

.111 

-.9 

AB— ii 

0.032 

0.303 

28 . 01 ; 

0.3768 

0.033 

O.03I; 

-2.9 

AB -3 

.032 

.388 

28 . 01 ; 

.7211 

.063 

.066 

-1.5 

AB-h 


0.313 


1.803 

0 . 0 U 60 

0 . 0 l ;22 

i -8.3 

AB -3 

■H 

.390 


1.082 

.Ol;la 

. 0 U 02 

: -8.9 
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TABLE II 


COMPARISON OF TESTED AND CALCDLATED DEFLECTIONS OF A SIMPLT 
SUPPORTED SQUARE SANDWICH PLATE, UNDER UNIFORM LOAD 

[^irminum faces with E = 10 x 10^ psi and cellular 
cellulose-acetate core ?d.th = 11,860 pal 
(reference 3); average absolute error, h‘9 percen'0 


Plate 

Thickness 
of faces 
(in.) 

Thickness 
of core 
(in.) 

Span 

(in.) 

Unit load 
on plate 
(psi) 

Test 

deflection 
(reference 3) 
(in.) 

Calculated 

deflections 

(in.) 

Percent 

error 


0.032 

0 .i ;93 

i 4 li .20 

m 

0.197 

0.186 

- 5.6 


.032 

.371 

Wi..20 


.316 

.306 

- 3.2 

AC-li 


O.liS? 

38 . Oil 

0.^768 

0.132 

0 . 121 ; 

"• 6.0 

AC -3 

■lEB 

.367 

38 . 0 i^ 

.5768 

.210 

.203 

-3.3 


0.032 

O.J189 

32. OU 

mmsm 

0.090 

0.083 

-7.8 

mi 

,032 

.36U 

32 . 0 li 

■^9 

. 1 U 6 

.137 

-6.2 

AC-li 

0.032 

O.U93 

mi 


0.073 

0.072 

-i.il 

AC -3 

.032 

.370 

HI 


. 08 l| 

.082 

-2.1; 


0.032 

0 .U 91 ; 


1.082 

0.0378 

0.0339 

-10.3. 


.032 

.367 


.7211 

. 0351 ; 

.0359 

l.h 


H 
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TABLE ni 

C0MPAEI50N OF TESTED AND CALCULATED DEFLECTIONS OF A SDIPLT 

SUPPORTED SQUARE SANDWICH ELATE UNDER UNIFORM LOAD 

[ at i-rnil mini faces with E - 10 X 10^ psi and corrugated 
paper honeycomb core with Qg = $8ij.0 psi (refer- 
ence 3); average absolute error, 6.^ percent] 


Plate 

Thickness 
of faces 
(in.) 

Thickness 
of core 
(in.) 

Span 

(in.) 

Unit load 
on plate 
(psi) 

Test 

deflection 
(reference 3) 
(in.) 

Calculated 

deflection 

(in.) 

Percent 

error 


0.032 

0.75^ 

Jlli.20 

0.7211 

0.110 

o.uh 

3.6 

AF-3 

.032 

.637 

. 14;. 20 

.7211 

.158 

.155 

-1.9 

AF-li 

0.032 

0.751 

38. OU 

0.7211 

0.072 

0.0653 

-9.3 

AF-3 

.032 

.632 

38, OU 

.7211 

.087 

.0890 

2.3 

AF-lt 

0.032 

0.752 

32.04 

0.7211 

0.041 

0.035 

-15 

aF -3 

.032 

.632 

32.04 

1 

.7211 

.054 

.047 

-13 

AF-k 

0.032 

0.7^3 

28.04 ' 

0.7211 

0.024 

0.021 

-12 

AP-3 

.032 

.633 

28.04 

.7211 

.030 

.029 

-3.3. 

AF-1; 

0.032 

0.75U 


2.884 

0.0391 

0.0384 

-1.8 

AF-3 

.032 

.636 


2.163 

.0368 

.0359 

-2.4 
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Figure 7-- Contour for S*(l3,0) - S*(a,0). 
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Figure 8.- Schematic representation for summing S(a,0) . 
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